What is claimed is: 



1. A two-stage method for characterizing a spatial arrangement 
among data points for each of a plurality of three--dimensional 
time series distributions comprising a sparse number of said data 
points, said method comprising the steps of: 

creating a first virtual volume containing a first three- 
dimensional time series distribution of said data 
points to be characterized; 



subdividing said first virtual volume into a plurality k of 
three-dimensional volumes, each of said plurality k of 
three-dimensional volumes having the same shape and 
size; 

providing a first stage characterization of said spatial 

arrangement of said first three-dimensional time series 
distribution of said data points comprising the steps 
of; 

determining a statistically expected proportion 0 of 
said plurality k of three-dimensional volumes 
containing at least one of said data points for a 
random distribution of said data points such that 
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Jc * © is a statistically expected number M of 
said plurality k of three-dimensional volumes 
which contain at least one of said data points if 
said first three-dimensional time series 
distribution is characterized as random; 



counting a number m of said plurality k of three- 
dimensional volumes which actually contain at 
least one of said data points in said first three- 
dimensional time series distribution;. 

statistically determining an upper random boundary 
greater than M and a lower random barrier less 
than M such that if said number m is between said 
upper random barrier and said lower random barrier 
then said first time series distribution is 
characterized as random in structure during said 
first stage characterization; 

providing a second stage characterization of said first 

three-dimensional time series distribution of said data 
points comprising the steps of; 
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when 0 is less than a pre-selected value, then 

utilizing a Poisson distribution to determine a 
first mean of said data points; 

when 0 is greater than said pre-selected value, then 

utilizing a binomial distribution to determine 
a second mean of said data points; 

computing a probability p from said first mean or from 
said second mean depending on whether 0 is 
greater than or less than said pre-selected value; 

determining a false alarm probability a based on a 
total number of said plurality Jc of three- 
dimensional volumes for said first three- 
dimensional time series distribution of said data 
points to be characterized; 

comparing p with a to determine whether to 

characterize said sparse data as noise or signal 
during said second stage characterization; and 

comparing said first stage characterization of said first 

three-dimensional time series distribution of said data 
points with said second stage characterization of said 
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first three-dimensional time series distribution of 
said data points to determine presence of randomness in 
said time series distributions. 

2. The method of claim 1, wherein if said first stage 
characterization of said first three-dimensional time series 
distribution of said data points indicates a random distribution 
and said second stage characterization of said first three- 
dimensional time series distribution of said data points 
indicates a signal, then continuing to process said data points. 

3. The method of claim 1, wherein if said first stage 
characterization of said first three-dimensional time series 
distribution of said data points indicates a random distribution 
and said second stage characterization of said first three- 
dimensional time series distribution of said data points 
indicates a random distribution, then labeling said' first three- 
dimensional time series distribution of said data points as 
random . 

4. The method of claim 1, further comprising utilizing the 
method steps of claim 1 for characterizing each of said plurality 
of three-dimensional time series distribution of said data 
points. 
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5. The method of claim 1, wherein said first three-dimensional 
time series distribution of said data points comprises less than 
about twenty-five (25) data points. 

6. The method of claim 1, wherein said upper random boundary- 
greater than M and said lower random barrier less than M are 
computed utilizing binomial probabilities. 

7. The method of claim 1, further comprising obtaining each of 
said plurality of three-dimensional time series distributions 
comprising said sparse number of said data points from a sonar 
system. 

8. The method of claim 1, further comprising obtaining each of 
said plurality of three-dimensional time series distributions 
comprising said sparse number of said data points from a radar 
system. 

9. The method of claim 1, further comprising determining said 
false alarm probability a based on a total number of said 
plurality k of three-dimensional volumes for said first three- 
dimensional time series distribution of said data points to be 
characterized wherein: 
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a 
a 



= 0.0U/A:>25, and 
= 0.05?/ it < 25. 



10. The method of claim 1, wherein. said step of comparing p with 

a to determine whether to characterize said sparse data as noise 

or signal during said first stage characterization is 

mathematically stated as: 

if p> a ^ NOISE, and 
if p<a=> SIGNAL. 



11. The method of claim 1, wherein said pre-selected value is 

equal to 0.10 such that if 

0 <0.10, then said Poisson distribution is utilized, and if 
0 >0.10, then said binomial distribution is utilized. 



12. The method of claim 1, wherein a total number Y of said data 

K 

points is given by Y = ^kNf^ , where 
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13. The method of claim 12, wherein said step of computing said 
probability p from said first mean further comprises utilizing 
the following equation: 

.1 



p = pL\<Zj=l—f= jexp(-.5x')dx 



where Z „ = ^ , 

where Y is said total number of data points, 
where, N is a sample size of said data points for each of a 
plurality of three-dimensional time series distributions, and 



//o="^T said first mean. 

k=0 



14. A method according to claim 13, wherein said step of 
computing said probability p from said second mean further 
conprises utilizing the following equation: 



/ X 1 

p = P(|za|<Z)=l — j= \exp{-.5x^)dx 

V2;r Hz,| 
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where = 



m±c-kO 



where c is a correction factor. 



15. The of claim 1, wherein k is determined from the relation 



kjj if Kj< Kfj , where 
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At is time interval .for collecting each of said plurality 
of three-dimensional time series distributions. 

Ay = inax(y)-min(y) where y is a magnitude of a first measure 
of said data points between a maximum and minimum value, and a 
second measure referred to as Z with magnitude 

AZ = max(Z)-min(Z) where Z is a magnitude of a second measure 
of said data points between a maximum and minimum value, and 

int is the integer operator. 
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